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Abstract 

£h ' 

Q \ The ac response of a slab of material with electrodynamic charac- 

^ ■ teristics E ~ j K , n > 0, is studied numerically. From the solutions 

^h^i of the nonlinear diffusion equation, the fundamental and higher-order 

components of the harmonic susceptibility are obtained. A large por- 
tion of the data for every k can be scaled by a single parameter, £ 
_£i/(/t+2) . ^«/( K + 2 ) jY)^ wriere i i s th e p er iod of the ac field at the 
surface, Hq is its amplitude and D is the slab thickness. This is, how- 
ever, only an approximate scaling property: The field penetration into 
a nonlinear medium is a more complex phenomenon than in the linear 
case. In particular, the susceptibility values are not uniquely defined 
by a set of only two parameters, such as k and £, while one parameter, 
i.e. t 1 / 2 /D, is sufficient to describe the electrodynamic response of a 
linear medium. 
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1. Introduction: Nonlinear diffusion and ac susceptibility of 
superconductors. 

The problem of nonlinear diffusion has recently attracted considerable 
attention in diverse fields of science. One example, which is considered in the 
present work, deals with the magnetization process of superconductors. In 
the case when the response of a superconductor to an applied ac field, H = 
H ■ exp(iut), is linear, the electrodynamic properties in the superconducting 
state can be described in terms of the complex conductivity, a = oy + i<Ji. 
In the flux-flow regime of high-T c superconductors, which occurs in a broad 
H — T range, it is usually justified to neglect the imaginary component of the 
complex conductivity. In that case, Maxwell's equations, VxB = 47r/c-j 
and cV x E = —dH/dt, provide the following result for the field penetration 
into a material filling half- space, x > 0: B = H ■ rj(x) ■ exp(iut), with 
rj(x) = exp(-Xx), and A 2 = 2i/5 2 , where 5 2 = c 2 / (27rua r ). For a thin plate 
of thickness 2D, with the ac field parallel to the large surface of the plate, 
the following equations hold [1]: 



. 1 sinhia) + sinia) „ 1 sin(a) — sinhia) 

Attx = -1 + rry- r^A^X = rpr- pr, (1) 

a cosn,{a) + cos{a) a cosh{a) + cos(a) 

where a = 2D/5. A maximum in x" results when 8(a r ,ou) is compa- 
rable to the sample size. Experimental results rarely show susceptibility 
curves which correspond to the flux-flow result of ohmic- like behaviour. 
Rather, the effects are most often nonlinear. The dependence of the mea- 
sured susceptibility curves on the excitation current and higher-order com- 
ponents in the harmonic susceptibilities are observed. In the limiting case 
of very strongly nonlinear response, the critical-state model may be used 
for the calculation of the ac susceptibility. Then, only one parameter is 
needed to construct the hysteresis curve for the magnetization, the field of 
the first full penetration to the sample center, H*. It is assumed that the 
harmonic susceptibility components, x' m an< ^ xl > are defined as Fourier 
components of the time-dependant magnetic hysteresis curve; M(t)/D = 
J2m (x'm ' cos(m ■ out) + x'L ' sin(m ■ cut)), where m is an integer. When, for 
instance, H < H* , Aitx'i = —(1 — H Q /2H*) is obtained then, 4:7ix' m — for 
every odd m > 1, and ^x'L = 2H /3nmH* for all odd m > 0. To deal with 
situations which are more relevant to the description of real experimental re- 
sults, it is necessary to investigate a nonlinear theory of the magnetic response 
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which would bridge the two limiting cases observed: the linear-response and 
the critical-state one. A fruitful approach to this problem is based on studies 
of the electrodynamic response of a medium characterized by a power-law 
current- volt age dependence, j = a(E) ■ E = a ■ E ■ (E / Eo) l ^ K+1 ^ , where 
k > 0. Using Maxwell's equations, the nonlinear diffusion equation describ- 
ing the penetration of fields into a slab of thickness 2D lying in the yz plane 
[1-3] can be derived as 



dfi__d_(dp_ 
dt dx \ dx 



8(3 



dx 



,t = -,x = —, (2) 

T X 



where (3 = B/Eq, xq = c/(47ro"o) and tq = l/(47ra"o). Recently, stud- 
ies of solutions of eq. (2) have been carried out by many authors [2-4]. 
The exact analytical description of the response of a superconductor to an 
abrupt change of external field has been given in references [1] and [2] and is 
compared with the results of non-logarithmic magnetization relaxation mea- 
surements on high-T c materials [5]. Various aspects of the ac response of 
superconductors has been studied as well by Dorogovtsev [6] and van der 
Beek et al. [7]. Recent results of Gilchrist and Dombre [8] can be compared 
with numerical results described in the present work. The distinctive feature 
of solutions of eq. (2) is that the flux-profile penetration resembles that in 
the models of the critical-state; When a field change is applied, the profile of 
perturbation spreads out from the surface towards the sample center but a 
region exists where the field distribution is unchanged inside. If the response 
to a field change of H is considered, the time after the front of the field 
change arrives to the center, t*, is given by, 

1 / r(l//t + l)r(3/2)\ K /D\ K+2 



r 2(k+1)(k + 2)\ r(l//c + 3/2) J \x J 
The initial magnetization, at t < t*, is given by: 

t \ 1 /(«+ 2 ) r(l//c + 3/2) \ 



- • (3) 



4ttM = -H ■ 1 - - • , , — ' ' : • (4) 

V r(i/« + 2) -r(i/2); v ; 

Equations (3) and (4) imply that a single parameter, £ = t 1 /( K+2 ) • 
can parametrize the short-time magnetization relaxation. It 
is informative to determine the extent to which this scaling relation is valid 
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with respect to the ac susceptibility (with the replacement of t and Hq by 
the ac field period and the field amplitude, respectively). 

2. Numerical modelling of nonlinear diffusion. 

Most of the calculations of the nonlinear diffusion process presented here 
have been performed on an array of dimension 50 x 200, containing mag- 
netic induction values, B, at 50 time intervals and 200 space-intervals [9]. 
The magnetic field at the surfaces of the sample, H Q ■ sin{uj ■ t), determines 
the boundary conditions. An average magnetic field < B(t) > in the sam- 
ple has been computed from the magnetic field distribution, every 50 time 
steps. Next, a Fourier time- analysis of < B(t) > has been performed and 
the coefficients of the fundamental and higher- order terms of the harmonic 
content have been found, 4:7ix' m =< B' m > and =< B'm except 

for the real component of the fundamental susceptibility, which is given by 
4t'i = — 1+ < B[ >. The method of computation and its results have been 
carefully tested. First, the magnetization relaxation process after an abrupt 
change of the external field has been simulated and numerical results were 
compared with the known exact analytical epressions derived by Koziol and 
de Chatel [2]. Then, the validity of the modelling of the ac-response in the 
limit of linear diffusion on the ac susceptibility, as given by eq. (1), was 
checked. It was confirmed also that the ac susceptibility converges towards 
the critical-state results for large k. The time range for which the response 
to the ac field becomes periodic (the initial response at short time does not 
satisfies this condition) was also investigated. In most cases it is safe to 
analyze the data taken after the initial 50000 steps in time evolution (this 
time depends on k and Hq). An additional, more reliable, criterion of stable 
periodicity is based on the criterion that the dc or second harmonic compo- 
nents are not found. At large values of H , the calculations become unstable 
abruptly. It is possible to overcome this difficulty but at significant expense 
in computation time (the computation of one susceptibility point requires an 
average of about 3 hours on an IBM-PC computer 486DX2-33MHz). There- 
fore, we have concentrated on performing calulations for a larger number of 
points at lower fields. 

3. Results and discussion 

The penetration of an alternating field resembles, in some ways, the 
response to an abrupt change of external field; the amplitude of field changes 
diminishes gradually in the material and, if the field amplitude at the surface 
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is not too large, there is no penetration to a volume separated by a certain 
distance from the sample surface. Whether the front of the flux profile in ac 
penetration propagates towards the center or not, is not an easy question to 
answer, since the initial very slow propagation which is observed might only 
be due to unstable initial conditions. Within the accuracy of calculations, 
the flux profile has a self-replicating shape of diminishing amplitude, with 
perfect periodicity in time at every point in space but with a phase shift 
which changes with the distance from the surface. The profiles obtained for 
one value of an ac field amplitude coincide with the profiles computed for 
another ac field amplitude, if the phase lag and spatial coordinates are shifted 
properly. Plots of x" versus x' shown in Figures 1 and 2 for different values 
of the nonlinearity parameter k converge to the limit of linear diffusion for 
k — > and to the limit given by the critical-state model for large k » 1. 




Figure 1. The x" versus x' plots of the ac susceptibility for different 
values of the nonlinearity parameter k and different periods of the ac field 
t, (n,t): (0.667, 25000, o ), (0.667, 5000, A), (2, 12500, • ), (2, 5000, □ ), 
(3, 10000, o ), (12, 6250, □ ). Solid lines represent the critical-state and the 
linear-response limits for a thin plate. The difference between the data for 
k = 2 obtained for two different frequencies of the ac field should be noted. 



5 



. 02 
.00 

CO 

X 

-0 . 02 
-0 . 04 



""-(7.05 -0.03 -0.01 0.01 

x 3 ' 

Figure 2. The third harmonic susceptibility x'l versus x'z compared 
to the critical-state result represented by the solid line, for the following 
nonlinearity parameter k and periods of the ac field t, (n,t): (0.667, 25000, 
• ), (0.667, 5000, O ), (1, 5000, D) and (2, 12500, o). 

An important feature of the present results is seen in the Figure 1; Sus- 
ceptibility points computed for different frequencies of the ac field but the 
same value of k, do not fall on the same curve. This is different from what 
it might be expected and seems to have been unnoticed in previous work 
[8]. In Figure 3, we show that a simple scaling of the susceptibility with the 
amplitude of the ac field holds for the data obtained in the range of incom- 
plete flux penetration, x ~ Hq^ k+2 \ In Figure 4, the real component of the 
first harmonic susceptibility is drawn as a function of H* I{k+2) ■ t^+^/D, 
for different values of the field amplitude H , period of the field and for a 
few sample sizes. This latter scaling method is not perfect; small differences 
in the slopes of the data computed for various frequencies is found. This 
effect may be explained by the fact that flux profiles have a shape which 
depends on the time of field penetration. One should expect that the param- 
eter £ = t 1 /(«+ 2 ) . jj*/( k + 2 ) i w jrj become an exact scaling variable only for 
the cases when all the parameters, t, Hq and D, are simultaneously scaled 
by a constant A in the following way: D — >• A • D, Ho — >• (A • iJ ) K /( K + 2 ) and 
i-.(A-t)V(«+2). 
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Figure 3. The real component of the third harmonic susceptibility as a 
function of the amplitude of the ac Geld, computed for an ac field of period 
equal to 25000. The slope of solid lines decreases for increasing values of 
n,which are the following: 0.667, 0.8, 1.2, 1.333, 1.667, 1.8, 2.333, 2.7, 3.35, 
3. 7. Similar scaling property is observed for the imaginary component of the 
third harmonic susceptibility and for higher-order components as well. 

4. Conclusions 

When the ac magnetic field does not penetrate to the sample center, the 
magnetic susceptibility is well described by a simple scaling relation: x ~ 



with £ = t 1 /^ 2 ) -Hq^^/D. For k values close to 0, the overall dependence 
of x"(x') closely resembles the dependence observed in the linear case. The 
conductivity, however, computed from x"(x') data by using the assumption 
that the linear theory holds, will lead to false information and yield strongly 
overestimated values. An experimental criterion for detecting nonlinearity 
would be the observation of the amplitude dependence or the existence of 
higher order harmonics in the ac response. The susceptibility values of a 
nonlinear medium are not uniquely defined by a set of two parameters only, 
such as k and £. For experimental purposes, however, treating £ as a scaling 
variable offers a sufficiently accurate method of testing for nonlinear proper- 
ties of materials. 
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Figure 4. Scaling of the flux penetration (4nx'i + 1) by a function 

gives the field amplitude, t gives the pe- 
riod and D gives the sample thickness. Each of the solid lines passes through 
the data corresponding to the following values of the nonlinearity parameter 
k: 0.667, 1, 2, 3 and 5, for lines with the smallest to largest slope. D is 
equal to 20 or 100 (there is no distinction between the symbols of the data 
corresponding to different values of D), while t is 2500 (o), 5000 (O), 6250 
(O), 10000 (A), 12500 (V), 25000 (•) and 100000 (D). 
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